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LIQUID MANAGEMENT IN LOW-GRAVITY USING BAFFLED ROTATING CONTAINERS 

INTRODUCTION 


Long duration space missions can require the location and control of large masses of liquids for long 

times. An extreme example of this problem is a proposed relativity experiment (GP-B * ) for which over 300 
kg of liquid helium will be vented during a period of one year, during which its configuration must be 

symmetric enough to keep acceleration effects at the design center of mass below 9.8 x 10"*® m/s^. It has 
been proposed that the helium, to be contained in a cylindrical annulus, be controlled by rotating the 
annulus, and providing baffles to help locate the helium. An ideal result is shown in Figure 1. The purpose 
of this paper is to examine the realizability of the configuration shown in Figure 1 . 

Without baffles the helium gas is likely to form a single bubble, which may arrange itself around the 
axis, or in an offaxis location. (For a full discussion of the unbaffled problem, see Reference 1.) From this 
it follows that there is likely to be a maximum permissable baffle spacing which will depend on the relative 
magnitude of the surface tension and centrifugal forces. As will be shown below, that proves to be the case; 
there is a baffle spacing beyond which the fluid cannot form the interface shown in Figure 1. The two 
possibilities are shown schematically in Figure 2. 

A second question is whether the interfaces will all remain at the same distance from the rotation 
axis. The layers of liquid are to be interconnected at the outside, so that withdrawal of fluid can be effected 
evenly from all the layers. Thus if there are pressure differences across the baffles, fluid will be forced from 
the high pressure layer to the low pressure layer. 

The pressure in the liquid is the difference between a centrifugal term and a surface tension term. 
If the interface moves out, both terms decrease. If the interface moves in, both terms increase. One cannot 
tell, a priori, whether the difference will increase or decrease. If the difference increases when the interface 
moves out of the configuration will be unstable because liquid will flow from the “shallow” layer to the 
“deep” layer. Figure 3 illustrates the situation for two layers. It will be shown below that this instability 
is possible at low rotation rates, and that it can arise for the presently-designed GP-B experiment. 

A third difficulty is the possible fluid dynamic instability of the interface itself. This problem is 
beyond the scope of this paper. 

Experimental analysis of possible fluid dynamical instability, and verification of the theoretical 

results to be given below, is to be desired. Unfortunately gravity levels of 9.8 x 10"*® m/s^ are not easily 
attainable. Thus the theory will be generalized to include an axial gravity field of arbitrary magnitude. A 
possible moderate gravity experiment will be discussed, and the results of a simple laboratory experiment 
will be given. Good agreement is demonstrated between calculation and observation for the latter. 



THE EQUATION OF THE INTERFACE 


Consider a liquid of density p ^ partially filling a cylindrical container of radius and length 

2(n+l)L. Let the cylinder contain n baffles, evenly spaced at intervals 2L (and neglect the baffle thickness.) 
Imagine the container to be rotating at co about its symmetry axis, which is antiparallel to a gravity field of 
strength g. 

Let the remainder of the container be filled with a gas of density pq < p^, and denote the inter- 
facial tension by y. Neglect compressibility. 

Construct a cylindrical r, <p, z coordinate system with r measured from the cylinder axis and z from 
the bottom of the cylinder. An interface shape 


r = rj(z) . (1) 

is to be found, compatible with a state of no relative motion — both gas and liquid are to be in a state of 
solid corotation with the container. Thus the pressure fields in the gas and liquid must be given by 

P G = p OG + 1 / 2 PG 0 * 2 * 2 ~ PG* 2 

( 2 ) 

P L = p OL + X < 2 PL^ 2 * 2 ~ Pl£ z 


where Pqq and Pql are constants. 

The interface need not be continuous because of the baffles. Indeed its derivatives are not 
continuous, and the interface itself can be continuous only in the absence of gravity. There will be n+1 
individual interface equations, the ith of which is 


r i = R i ~ ^i( z > 


(3) 


where is a constant, representing the maximum distance from the rotation axis to the interface. 

Figure 4 shows a single layer, and introduces some additional notation: is the contact angle and 

z^- and z- P are the distances between the f ■ = 0 position and the upper and lower bounding baffles. 

The jump in pressure across each interface is equal to 7K', where K' denotes the mean curvature. 
Using standard formulas from differential geometry [2] it can be shown that 


K' = i4^ 

r r 3 dz 2 


(4) 


2 



where 


A = 



2 1/2 


(5) 


and the subscript i has been suppressed. 

The interface equation in the ith layer can be written in terms of local dimensionless coordinates 
Xj and y^, defined by 


z = 2(i-l)L + z i B + R i x i 
r=R i = R i [l-y i (x i )] 


( 6 ) 


The result is 


1 _ 

A 


1+ ^4 F i y i (2 - y i } 


- GjXj = 1 + Cj 


(7) 


where the prime denotes differentation with respect to argument, 
Fj = pco 2 R i 3 /7 Gj = pgRj 2 /7 


( 8 ) 


p = PL - Pq is the density contrast and 



(9) 


is the axial curvature at the point x^ = 0. 

Because the various liquid segments are connected at their outer edges, the liquid pressure must be 
continuous there. This reduces to the condition that 


r,{ 


n =— U+C i +l/2F i -G i 2(i 


• zBi 1 
*i R i J 


( 10 ) 


be independent of i. (In the absence of gravity, one solution is the intuitive one that Rj = R, a constant.) 
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Except for the condition (10), the subscript i is not necessary and will be suppressed below where 
possible. 

An interface shape is to be found by integrating equation (9) subject to the boundary condition that 


y' = cot0 at x = z^/R 

y' = - cot0 at x = -z®/R 


(ID 


and the compatibility condition that z^ - z® = 2(L/R). 

A rescaling of equation (7) assuming y -+• ey, z ->• ez, e 2 G « 1 , e e F = 0(1) leads to the usual equa- 
tion for a planar interface under a perpendicular gravity field [3], which can be integrated directly. That 
case is not relevant to the present inquiry, and it is necessary to integrate equation (7) numerically. That is 
most easily done by rewriting the equation as a pair of first order equations 


y' = Y 


Y' = (l+Y 2 ) 



[ 1 EC-Gz+Fy ( 1 - 1 /2y) ] 


1 

l^y 


( 12 ) 


and then applying a fourth order Runge-Kutta scheme [4] to this system. It is only necessary to note that, 
for small contact angles, Y becomes large near the wall. This necessitates a continual adjustment of the 
interval in z. The easiest way to do this is to fix the interval As along the interface and to write the interval 
in z, h, as 


h = As/(1 + Y 2 ) 1 / 2 


(13) 


The actual integration is performed, for a given F, by starting at z = 0 (where y = 0 = y') with a given value 
of C and integrating first in the forward direction until Y = cot0 and then, starting again at z = 0, integrating 
in the negative direction until Y = -cot0. Finally the compatibility conditions are applied, in effect deter- 
mining L/R a posteriori. 


RESULTS AND DISCUSSION: STABILITY, AN EXPERIMENT AND GP-B 


Interface profiles, and results derived from them, are obtained by choosing co and g, picking values 
of R and C, and then integrating equations (12), as described above, until the boundary condition (1 1) are 

met. The values of z^,z® and 2L are recorded, as well as the maximum values of y, the capillary rise at each 
end, y^ and y®. For fixed co, g and R (hence F and G z^, z®, y^ and y® all increase as C decreases. For any 
F and G there appears to exist a minimum value of C, at which the larger y^ or y® reaches unity, which 
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corresponds to a capillary rise reaching the rotation axis. Decreasing C below this value leads to an apparent 
change in sign in y" before y' can reach its boundary value. Because one feels that the interface slope should 
be monotonic, and because the maximum capillary rise cannot exceed the maximum radius R, this critical 
value of C, C m j n and the value of L, L max associated with it are considered to be real limits. Specifically 

2L max is assumed to be the maximum baffle spacing for which an interface can form. 

Figure 5 shows the dimensionless maximum baffle spacing D = 2(L/R) as a function of F for 6 = 
0 = G. Increasing the contact angle increases D. Figure 6 shows the effect of adding gravity, a decrease in D. 

The second type of instability is that outlined in the introduction, for which the dependence of 
pressure on radius is such that a perturbation of the equilibrium position tends to grow. If the (i+l)st slice 
is moved outward, so that Rj + j -*■ Rj + j + ARj + j, and the ith slice is moved inward, so that Rj -*■ Rj - ARj, 

then the pressure drop from i+1 to i is given by 



AR 


i+1 


+ 



AR; 


and if that is positive, then the displacement will grow and the pair is unstable. 

It is necessary to find the rate of change of the liquid pressure with R, holding co, g and L fixed. 
To that end, write the dimensional pressure 


P L = P 0L + < P 0G ~ P 0G> + !/2 PL w2 r2 ~ p jJ> z 


1+C; 


= PqG * T ~ ! / 2 P^ 2 Rj 2 + Pg 2(i-l) L + Zj B + 1/2 p L co 2 r 2 - P L g 2(i-l) L + Zj B + Xj 


= 1/2 p L w 2 r 2 - p L g 2(i-l) L + + xj + Pq G - yllj . 


B 


(14) 


The pressure in question is that at x = -Zj B or x = . In either case x^ + Z| B is independent of Rj and 


9P 

”aR 


L an: y f aq a Zi B 

| 1+c i- F i- R i9R- +G i«r ■ 


(15) 


In any such virtual displacement of the interfaces, volume must be conserved. Thus AR- cc R^‘ ^ , and 
the criterion for the interchange instability is 


1 9n i 


1 dn i+l 

^__ + _ <0 


Rj 3R R i+1 3R 


(16) 
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If G = 0 the situation is simplified dramatically because the equilibrium position has Rj = R] + ] , C] = C] + j, 
F- - Fj + i and the neutral curve is given by 


2 25L-1 

9R 


+ C-F-R 


9 C 
9R 


= 0 


(17) 


In the case of zero g, where equation (17) applies, C has been found to be well represented in the 

form 


C = A(F) (R/L) + B(F) 


(18) 


so that 


F = r +3 |- [A ' (R/L)+B,] 


(19) 


where ' denotes derivative with respect to argument, and equation ( 1 7) can be rewritten 

1 - F + B - 3F[A'(R/L) + B'] = 0 . (20) 


that curve is shown in Figure 7. One notes that this is the limiting stability curve for F < 1.2, and that it 
gives a stringent limit for F < 1 ,0. 

The situation for nonnegligible G is sufficiently complicated that a general analysis seems counter- 
productive. A possible application to moderate g (a low g profile term on an aircraft will be explained as 
an example. 

Consider a two layer model for which y = 0.022 N/m, p = 790 kg/m^, 2L = 20 mm, co = 1 s'* and 

g = 0.098 m/s^. Figure 8 shows II] and II 2 as a function of R, where 1 denotes the lower layer and the 2 
the upper. 

Because all the independent variables other than the void volume have been fixed, R is a function of 
gas volume only, and the position of any configuration is determined by that void volume. Denote that by V. 
For V sufficiently small, the system cannot form an interface at the given spacing, D. Let this critical volume 
be Vj. In this case Vj ~ 9 ml, corresponding to R ~ 10 mm. Note that the II, R curves do not extend to the 

left of this critical point. 

There is a second critical volume, V 2 , below which two interfaces are unstable to the interchange 
instability. Such a circumstance is shown by the configuration B, which has a total void volume, Vg, less 

than that of the configuration A,V^; x For the state B, 9II/9R < 0 and 9112/ ^ ^ the calculation 
shows that 
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1 an l 1 3n 2 

r7 ~3R~ + rJ “aR" 


so that the configuration is unstable. The only other configuration with the same volume is one for which 
the lower layer fills in to the center and the upper layer becomes an apparent first layer. This configuration 
is that of B', and the theory given in this paper predicts that the B' configuration is that realized. 

For volumes greater than that associated with the configuration A, both Bflj/BR and 9E[j/3R are 

positive. Such configurations will be stable. So will the alternative configuration A!, in which the lower 
layer is full and upper acts as a lower. The present theory cannot answer the question of the state attained 
as V goes from Vg to V^. 

One motivation for adding gravity was the hope that useful experiments could be done in the labora- 
tory. Ideally one would want small G and moderate F to attempt to model the proposed space application, 
for which G « 1 and F = 0(1). That proves to be impossible for reasonable values of R and co. If p, g and 
7 are fixed, then 


R= [G7/pg] 1/2 


and 


co = [F 7 /pR 3 ] */ 2 = (pg 3 / 7 ) 1/2 FG' 3 / 2 


and, to use ethanol as an exemplary fluid, if G = 0.1 and F = 1, then R = 537 pm and co = 1.8 x 10 3 s‘*. 

While a model experiment is impossible, it seemed useful to explore the basic premises of the model 
in the laboratory. To that end, a set of experiments in a shallow disk were performed. All that could be 
observed easily was the location of the top interface (when the liquid intersects the upper boundary) and 
whether the liquid intersects the lower boundary. (If 2L exceeds the maximum value for a given F, the 
liquid will intersect only the top surface.) 

The experiments were carried out in a small petri dish (diameter = 48.06 mm, depth = 8.59 mm) 
mounted, using double sided tape, to the top of a Genisco model C-181 turntable (maximum rotation rate 

22 s" 1 ). Centering accuracy was better than 0.8 mm. Turntable speed was accurate to better than 5 parts 
per thousand in the range of interest. The working fluid was ethanol (Pharmco 200 proof ethyl alcohol) 
colored with Higgens India Ink at 2 drops for 25 ml ethanol. Textbook [3] values for density (800 kg/m 3 ) 
and surface tension (0.022 N/m) were used. A contact angle of zero was assumed. 

The measured dependent variable was the apparent intersection radius of the fluid with the upper 
boundary as a function of rotation rate and air volume. The observation was made visually by watching the 
position of the intersection against a marked upper surface, marked with concentric circles spaced at 2 mm 
intervals. Errors arise from nonconcentricity and parallax (because the circles are on the upper surface of the 
cover, to avoid imposing surface drag on the contact line.) These errors are estimated to be 0.5 mm. The 
mean position can be estimated somewhat more precisely than 0.5 mm. Positions shown are to 0.1 mm, 
which may be optimistic. 
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There is an additional systematic error which may lead one to overestimate the radius. This arises 
because the visible intersection line is not the true intersection line, but is some point at which the coloring 
is dark enough to see. The nearer the contact angle is to zero, the more pronounced this effect will be. 

The most difficult parameter to control in the experiment is the air volume. Several methods were 
tried. The final procedure was to fill the container as nearly completely as possible, the last few ml being 
added through a small hole in the center of the cover through a syringe while the container was spinning 

at 22 s' 1 , to center the air bubble; ( 2 ) withdraw the desired amount of fluid with the container stationary; 
(3) plug the hole with paraffin. The last step is necessary because, when the hole is unplugged, the ethanol 
is sufficiently volatile to evaporate at about 1.2 ml/hr (measured with the container stationary). 

Results are shown in Figure 9 for air volumes of 2, 4, and 6 ml. The symbols denote data and the 
solid lines joining open circles are calculated results. The calculation predicts-bottom exposure for the 6 ml 

case at rotation rates above 18 s' 1 , and the observations are consistent with the prediction. The four right- 
most symbols on the upper curve showed a clearly exposed bottom. The next two were ambiguous. All 
the others, on all the curves, showed the bottom covered. 

Finally, the consequences of this analysis for the proposed GP-B experiment mentioned in the 
introduction should be considered. The container for that application is a cylindrical annulus with an inner 
radius 0.18 m and R„ = 0.54 m. The overall length is 2.94 m. Reasonable values for the physical constants 
c 1 

of liquid helium are [1] : p = 145 kg/nr and 7 = 5.3 x 10 ' 1 N/m. The design rotation rate is to be 0.01 

s' 1 . The dewar is to be initially 90 percent full of liquid helium, This provides an estimate of R > 0.242 m, 
and gives values of G and G, 


F « 0.388 


G ss 1.56 x 10' 5 > 


the latter based on the necessity of keeping local gravity to 10' 111 x normal gravity. 

For these values of F and G one can find the maximum baffle spacing for which the helium does 
not reach the surface of the inner annulus. That result is 2L = 155 mm. This, however, is not the limiting 
value, because the interchange instability is dominant at this small value of F. By repeated evaluation of 
equation (1 7) at varying values of 2L, the critical values is found to be 2L = 18.5 mm. 
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Figure 6. Neutral stability curves for nonzero gravity: D,F pairs above each curve cannot form interfaces at G 
the value shown on the curves. Points denote numerical calculations. 
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